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Abstract 

The Gauss - Bonnet invariant is one of the most promising candidates for 
a quadratic curvature correction to the Einstein action in expansions of su- 
persymmetric string theory. We study the evaporation of such Schwarzschild 
- Gauss - Bonnet black holes which could be formed at future colliders if 
the Planck scale is of order a TeV, as predicted by some modern brane 
world models. We show that, beyond the dimensionality of space, the cor- 
responding coupling constant could be measured by the LHC. This opens 
new windows for physics investigation in spite of the possible screening of 
microphysics due to the event horizon. 



PACS Numbers: 04.70.Dy (Quantum aspects of black holes), ll.25.-w (Strings 
and branes), 13.90. +i (phenomenology of elementary particles) 



1 Introduction 

It has recently been pointed out that black holes could be formed at future col- 
liders if the Planck scale is of order a TeV, as is the case in some extra-dimension 
scenarios PQEj- This idea has driven a considerable amount of interest (see e.g. 
jSj). The same phenomenon could also occur due to ultrahigh energy neutrino 
interactions in the atmosphere jlj. Most works consider that those black holes 
could be described by the D-dimensional (D > 5) generalized Schwarzschild or 
Kerr metrics The aim of this paper is to study the experimental consequences 
of the existence of the Gauss - Bonnet term (as a step toward quantum gravity) if 
it is included in the D-dimensional action. This approach should be more general 
and relies on a real expansion of supersymmetric string theory. In Section 2, the 
basics of black hole formation at colliders and the related cross sections are re- 
minded. The details of the mult i- dimensional Gauss - Bonnet black hole solutions 
and their thermodynamical properties are given in Section 3. The flux compu- 
tation and the main analytical formulae are explained in Section 4. It is shown 
in Section 5 that the Gauss - Bonnet (string) coupling constant can be measured 
in most cases, together with the dimensionality of space. Finally, some possi- 
ble consequences and developments, especially with an additional cosmological 
constant, are discussed. 

2 Black hole formation at colliders 

The "large extra dimensions" scenario |0] is a very exciting way to address ge- 
ometrically the hierarchy problem (among others), allowing only the gravity to 
propagate in the bulk. The Gauss law relates the Planck scale of the effective 
4D low-energy theory Mpi with the fundamental Planck scale through the 
volume of the compactified dimensions, Vd-a, via: 



It is thus possible to set Mo ~ TeV without being in contradiction with any 
currently available experimental data. This translates into radii values between 
a fraction of a millimeter and a few Fermi for the compactification radius of the 
extra dimensions (assumed to be of same size and flat, i.e. of toroidal shape). 
Furthermore, such a small value for the Planck energy can be naturally expected 
to minimize the difference between the weak and Planck scales, as motivated by 
the construction of this approach. In such a scenario, at sub-weak energies, the 
Standard Model (SM) fields must be localized to a 4-dimensional manifold of 
weak scale "thickness" in the extra dimensions. As shown in [6:, as an example 
based on a dynamical assumption with D=6, it is possible to build such a SM 
field localization. This is however the non-trivial task of those models. 

Another important way for realizing TeV scale gravity arises from properties 
of warped extra-dimensional geometries used in Randall- Sundrum scenarios jjj. 
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If the warp factor is small in the vicinity of the standard model brane, particle 
masses can take TeV values, thereby giving rise to a large hierarchy between 
the TeV and conventional Planck scales [2J IE]- Strong gravitational effects are 
therefore also expected in high energy scattering processes on the brane. 

In those frameworks, black holes could be formed by the Large Hadron Col- 
lider (LHC). Two partons with a center-of-mass energy ^/s moving in opposite 
directions with an impact parameter less than the horizon radius r + should form a 
black hole of mass M ~ \fs with a cross section expected to be of order a ~ nr 2 ,. 
Thoses values are in fact approximations as the black hole mass will be only a 
fraction of the center-of-mass energy whose exact value depends on the dimen- 
sionality of the spacetime and the angular momentum of the produced black hole 
[HI [TUJ . Furthermore, suppression effects in the cross section should be consid- 
ered and are taken into account in the section 5 of this paper. Although the 
accurate values are not yet known, a semiclassical analysis of quantum black hole 
formation is now being constructed and the existence of a closed trapped surface 
in the collision geometry of relativistic particles is demonstrated. To compute 
the real probability to form black holes at the LHC, it is necessary to take into 
account that only a fraction of the total center-of-mass energy is carried out by 
each parton and to convolve the previous estimate with the parton luminosity 
P]. Many clear experimental signatures are expected in particular very high 
multiplicity events with a large fraction of the beam energy converted into trans- 
verse energy with a growing cross section. Depending on the value of the Planck 
scale, up to approximately a billion black holes could be produced at the LHC. 



3 Schwarzschild - Gauss - Bonnet black holes 



The classical Einstein theory can be considered as the weak field and low energy 
limit of a some quantum gravity model which is not yet built. The curvature 
expansion of string gravity therefore provides an interesting step in the modelling 
of a quasiclassical approximation of quantum gravity. As pointed out in [TTj . 
among higher order curvature corrections to the general relativity action, the 
quadratic term is especially important as it is the leading one and as it can affect 
the graviton excitation spectrum near flat space. If, like the string itself, its slope 
expansion is to be ghost free, the quadratic term must be the Gauss - Bonnet 
combination : Lqb = R^ a pR^ aP ~ 4R a pR af3 + R 2 . Furthermore, this term 
is naturally generated in heterotic string theories and makes possible the 
localization of the graviton zero- mode on the brane . It has been successfully 
used in cosmology, especially to address the cosmological constant problem (see 
e.g. [TJ] and references therein) and in black hole physics, especially to address 
the endpoint of the Hawking evaporation problem (see e.g. ^5] and references 
therein). We consider here black holes described by such an action : 



S 
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where A is the Gauss - Bonnet coupling constant. The measurement of this A 
term would allow an important step forward in the understanding of the ultimate 
gravity theory. Following ^H] > we assume the metric to be of the following form : 

ds 2 = -e 2u dt 2 + e 2a dr 2 + r 2 h ij dx i dx j 

where v and a are functions of r only and hi^dx % dx^ represents the line element of 
a (D — 2)-dimensional hypersurface with constant curvature (D — 2)(D — 3). The 
substitution of this metric into the action [TT] leads to the following solutions : 

r 2 

e 2u = e -2a = 1 + 



2X(D -3)(D -4) 



i , t 32tt^ GX{D -3)(D -4)MT( 



D-l- 



(D - 2)r D -! 

The mass of the black hole can then be expressed fTJ QS] m terms of the horizon 
radius r + , 

(D-2)«Jr^( A( g -3)( C -4) 

8irGr (s=i) V r i 

where T stands for the Gamma function. The temperature is obtained by the 
usual requirement that no conical singularity appears at the horizon in the eu- 
clidean sector of the hole solution, 

rp J_ ( -2«y , (D-3)r 2 + + (D-5)(D-4)(D-3)\ 

BH 4n [e } W=r+ 4nr + (rl + 2X(D-4)(D-3)) ' 

In the case D = 5, those black holes have a singular behavior ^H] and, depending 
on the value of A, can become thermodynamically unstable or form stable relics. 
For D > 5, which is the only relevant hypothesis for this study (as D = 5 would 
alter the solar system dynamics if the Planck scale is expected to lie ~TeV), a 
quantitatively different evaporation scenario is expected. Figure H shows the ratio 
of the temperatures with and without the Gauss - Bonnet term for different values 
of D and A. It should be pointed out that the non-monotonic behavior makes 
an unambiguous measurement quite difficult and requires to take advantage of 
the full dynamics of the evaporation. The next sections focus on this point to 
investigate the A parameter reconstruction. 



4 Flux computation 

Using the high-energy limit of multi-dimensional grey-body factors J7j, the spec- 
trum per unit of time t and of energy Q can be written, for each degree of freedom, 
for particles of type i and spin s as: 

d 2 ^ _ ^ 2 (^)^m)riQ 2 

dQdt _ (_i)2 S 



3 




Figure 1: Ratio of the temperatures with and without the Gauss - Bonnet term 
for D = 6, 7, 8, 9, 10, 11 (from up to bottom in the low mass region) as a function 
of mass with A = 1 TeV" 2 (up) and A = 0.01 TeV~ 2 (down). 



This is an approximation as modifications might arise when the exact values of 
the greybody factors are taken into account due to their dependence, in the low 
energy regime, on both the dimensionality of the spacetime and on the spin of the 
emitted particle. Fortunately, as demonstrated in the 4-dimensional case ^H|, the 
pseudo-oscillating behaviour induces compensations that makes the differences 
probably quantitatively quite small. As shown in the previous section, as long 
as D > 5, the horizon radius r + cannot be explicitly given as a function of the 
mass and, to compute the experimental integral spectrum diVj/dQ, the following 
change of variable is convenient : 
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Nf and iV& being the total fermionic and bosonic degrees of freedom. The mean 
number of emitted particle can then be written as 



Ni 
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where ri n u+ is the initial horizon radius of a black hole with mass Mi n u and, 
interestingly, the ratio of a given species i to the total emission is given by : 

Ntot ~ \N f + N tot 

where a s is 1 for bosons and is 3/4 for fermions and ^ is the number of internal 
degrees of freedom for the considered particles. The mean number of particles 
emitted by a Schwarzschild - Gauss - Bonnet black hole ranges from 25 to 4.7 
depending on the values of A and D, for M D ~ 1 TeV and M init ~ 10 TeV. Those 
values are decreased to 5 and 1.05 if M init is set at 2 TeV. Figure El shows the flux 
for different values of A and D. Although some combinations seem to be strongly 
degenerated, the next section shows that in any case the values of A and D can 
be well reconstructed. 





Figure 2: Integrated flux as a function of the total energy of the emitted quanta for 
an initial black hole mass M = 10 TeV. Upper left : A = 0, D = 6, 7, 8, 9, 10, 11. 
Upper right : A = 0, 5 TeV~ 2 , D = 6,7,8,9,10,11. Lower left : D = 6, A = 
0.1, 0.5, 1, 5, 10 TeV" 2 . Lower right : D = 11, A = 0.1, 0.5, 1, 5, 10 TeV~ 2 . 
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5 String coupling constant measurement 



To investigate the LHC capability to reconstruct the fundamental parameter A, 
we have fixed the Planck scale at 1 TeV. Although a small excursion range around 
this value would not change dramatically our conclusions, it cannot be taken much 
above, due to the very fast decrease of the number of formed black holes with 
increasing Mb- Following [ X J , we consider the number of black holes produced 
between 1 TeV and 10 TeV with a bin width of 500 GeV (much larger than 
the energy resolution of the detector), rescaled with the value of r + modified by 
the Gauss - Bonnet term. For each black hole event, the emitted particles are 
randomly chosen by a Monte-Carlo simulation according to the spectra given in 
the previous section, weighted by the appropriate number of degrees of freedom. 
The Hawking radiation takes place predominantly in the S-wave channel so 
bulk modes can be neglected and the evaporation can be considered as occurring 
within the brane. As the intrinsic spectrum diVj/dQ is very strongly modified by 
fragmentation process, only the direct emission of electrons and photons above 
100 GeV is considered. We have checked with the Pythia [20] hadronization 
program that only a small fraction of directly emitted 7-rays and electrons fall 
within an hadronic jet, making them impossible to distinguish from the back- 
ground of decay products. Furthermore, the background from standard model 
Z(ee)+jets and 7+jets remains much lower than the expected signal. The value 
of the Planck scale is assumed to be known as a clear threshold effect should 
appear in the data and a negligible uncertainty is expected on this measurement. 
For each event, the initial mass of the black hole is also assumed to be known as 
it can be easily determined with the full spectrum of decay products (only 5% of 
missing energy is expected due to the small number of degrees of freedom of neu- 
trinos and gravitons). The energy resolution of the detector is taken into account 
and parametrized [2T| as a/E = ^Ja? /E + b 2 with a m 10%v / GeV and b « 0.5%. 
Unlike pQ, we also take into account the time evolution of the black holes and 
perform a full fit for each event. Once all the particles have been generated, 
spectra are reconstructed for all the mass bins and compared with theoretical 
computations. The values of D and A compatible with the simulated data are 
then investigated. Figure El shows the x 1 /d.o.f. for the reconstructed spectra 
for 2 different couples (A [ TeV~ 2 ], D) = (l,10) and (A [ TeV~ 2 }, £>) = (5,8). The 
statistical significance of this x 2 should be taken with care since a real statistical 
analysis would require a full Monte-Carlo simulation of the detector. Neverthe- 
less, the "input" values can clearly be extracted from the data. Furthermore, it 
is important to notice that for reasonable values of A (around the order of the 
quantum gravity scale, i.e. around a TeV -2 in our case) it can unambiguously be 
distinguished between the case with and the case without a Gauss - Bonnet term. 
Table n summarizes the LHC reconstruction capability requiring the x 2 /d.o.f. to 
remain smaller than 2x 2 nin / d.o.f . where Xmin/ 'd.o.f '. corresponds to the "physi- 
cal" case (i.e. A = Xi npu t and D = D input ). This is quite conservative and should 
translate into high confidence levels which would require a much more detailed 
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modelling of the detector to be accurately computed. For each set of parame- 
ters, the cross section has been taken as 7rr^_, 7rr+/10, 7rr^_/100 and 7rr^/1000 to 
account for uncertainties on the production process for D > 4 with a non-zero 
impact parameter. Based on the methods developed by Penrose and D'Eath & 
Payne and on the hoop conjecture ^U], several estimates have been derived 
and confirm the formation of an apparent horizon. The wide range investigated 
in this study should account for all physical cases. 




Figure 3: Upper part : values of the x 2 /d.o.f . for the reconstructed spectra as a 
function of D and A for "input" values A = 1 TeV~ 2 and D = 10 ; the right side 
shows rectangles proportional to the logarithm of the x 2 /d.o.f . Lower part (left 
and right) : values of the x 2 /d.o.f . for the reconstructed spectra as a function 
of D and A for "input" values A = 5 TeV~ 2 and D = 8 ; the right side shows 
rectangles proportional to the logarithm of the x 1 /d.o.f . 
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Table 1: Reconstructed values for D and A (TeV~ 2 ) as a function of the "real" 
input values requiring x 2 < 2Xmin- The first line assumes a = 7rr 2 , the second 
line o = 7rr 2 /10, the third line a = 7rr 2 /100 and the fourth line a = 7rr 2 /1000. 



6 Discussion 

In case the Planck scale lies in the TeV range due to extra dimensions, this 
study shows that, beyond the dimensionality of space, the next generation of 
colliders should be able to measure the coefficient of a possible Gauss - Bonnet 
term in the gravitational action. This would allow an important step forward 
in the construction of a full quantum theory of gravity. It is also interesting to 
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notice that this would be a nice example of the convergence between astrophysics 
and particle physics in the final understanding of black holes and gravity in the 
Planckian region. 

Nevertheless, those results could be improved and refined in several ways. 
First, the endpoint of the Hawking evaporation process is still an unsolved prob- 
lem. In this paper, we have considered that the time integral of the instantaneous 
spectrum is valid up to the total disappearance of the black hole. Although usu- 
ally a good approximation (as most particles are emitted at masses close to the 
initial mass), this can become a serious problem if the number of extra dimensions 
is high. In such cases, the mean number of emitted particles can be very small 
and even smaller than one. The spectrum therefore must be truncated properly. 
A possibility could be to add a Heavyside function to ensure energy conservation 
while keeping the same probability distribution, as suggested in [22], but a full 
understanding of the phenomenon would be required as the analytical formulae 
derived in this work would not stand anymore. 

Then, as studied in [1^1221, a cosmological constant could also be included 
in the action. On the theoretical side, this would be strongly motivated by the 
great deal of attention paid to the Anti-de Sitter and, recently, de Sitter / Con- 
formal Field Theory (AdS and dS /CFT) correspondences. On the experimental 
side, this would open an interesting window as there is no unambiguous relation 
between the D-dimensional and the 4-dimensional cosmological constants. 

Finally, it would be very interesting to extend this study to Kerr - Gauss - 
Bonnet black holes as the holes possibly produced at colliders are expected 
to be spinning. Although qualitatively equivalent, the results are expected to be 
quantitatively quite different and probably more realistic. 
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